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Abstract. The drag and diffusion coefficients of heavy quarks (HQs) have been 
evaluated in the pre-equilibrium phase of the evolving fireball produced in heavy ion 
collisions at RHIC and LHC energies. The KLN and classical Yang-Mills spectra have 
been used for describing the momentum distributions of the gluons produced just after 
the collisions but before they thermalize. The interaction of the HQs with these gluons 
has been treated within the framework of perturbative QCD. We have observed that the 
HQs are dragged almost equally by the kinetic.ally equilibrated and out-of-equilibrium 
gluonic systems. We have also noticed that the HQs diffusion in the pre-equilibrium 
gluonic phase is as fast as in the kinetically equilibrated gluons. Moreover, the diffusion 
is faster in the pre-equilibrium phase than in the chemically equilibrated quark-gluon 
plasma. These findings may have significant impact on the analysis of experimental 
results on the elliptic flow and the high momentum suppression of the open charm and 
beauty hadrons. 
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1. Introduction 

The primary intent of the ongoing nuclear collision programmes at the Relativistic 
Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) energies is to create a 
new state of matter called Quark Gluon Plasma (QGP). The bulk properties of this 
state is governed by the light quarks (q) and gluons (g). The heavy quarks [HQs = 
charm (c) and beauty (b)] are considered as efficient tools to probe the early state of 
the system (see [BEIE] for a review) mainly for the following reasons: (i) HQs are 
produced at very early stage; (ii) The probability of creation and annihilation of the 
HQs during the evolution of the fireball is small. Hence the HQs can witness the entire 
evolution of the system as they are created early and survive the full evolution without 
being annihilated and/or created. 

The data on the suppression of the charm quark at large momentum (Raa(pt)) 
|2 El El U\ and their elliptic flow (w 2 ) [3 E] have been analyzed by several authors (see 
e.g. [2 E.» Ej] and references therein) to characterize the system formed in HIC at RHIC 
and LHC collisions. 

It has been argued that the thermalization times of HQs are larger than the light 
quarks and gluons PHD] by a factor ~ M/T [lTj , where M is the mass of the HQs and 
T is the temperature of the medium. Therefore, the evolution of the HQs in the thermal 
bath of light quarks and gluons can be treated within the ambit of Brownian motion, 
although a detailed investigation of this problem within the framework of Boltzmann 
equation has revealed that the evolution of charm quarks as a Brownian particle requires 
some corrections [ 12 ] (see also psi ei 112 m ei ns]). 

Several theoretical attempts have been made to study the evolution of the HQs 
within the framework of Fokker Plank equation [2 [Li] [191 EOl E2 E2J El] E2 EH EH] 127] 
es EH m eh m eh However, the roles of pre-equilibration phase have been ignored 
in these works. An attempt has been made in the present work to study this role. 
This is important because the HQs are produced in the hard collisions of the colliding 
partons of the nuclei and they inevitably interact with the pre-equilibrium phase of 
the bulk matter. The effect of the pre-equilibrium phase might be more significant 
for low-energy nuclear collisions: for example, in the case Au+Au collisions at RHIC 
energy simulations show that equilibration is achieved approximately within 1 fm/c, 
while the lifetime of the QGP phase turns out to be about 5 fm/c [12 35]; hence the 
lifetime of the out-of-equilibrium phase is approximately 20% of the total lifetime of 
the QGP. The lifetime of the QGP phase is of about 10 fm/c and the duration of the 
pre-equilibrium phase is about 0.4 fm/c for LHC collision condition. Although, this 
may suggest the dwindling role of the pre-equilibrium phase for collisions at higher 
energies, estimates of drag and diffusion coefficients of HQs done in this work indicate 
non-negligible contributions of the pre-equilibrium phase. The motivation of this work 
is to estimate the drag and diffusion coefficients of the HQs interacting elastically with 
the non-equilibrated gluons constituting the medium. 

The paper is organized as follows. In the next section we discuss the formalism used 
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to evaluate the drag and diffusion coefficients of the heavy quarks in the preequilibrium 
stage. In Section III we summarize the out-of-equilibrium initial conditions we 
implement in the calculations. Section IV is devoted to present the results and section 
V contains summary and discussions. 


2. Formalism 


The Boltzmann Transport Equation (BTE) describing the evolution of the HQs in the 
pre-equilibrated gluonic system can be written as: 

~df 


d p d -> d 

- . -|_ F ■ — 

dt E dx dp 


dt 


( 1 ) 


f(x,p,t) = 

/ it 

collisions 

For binary interaction the collisional integral appearing in the right hand side of BTE 
can be written as: 

~df 


dt 


= / d 3 k[w(p + k, k)f(p + k) — w(p, k)f(p)\. 


( 2 ) 


collisions 


where w(p, k) is the rate of collision which encodes the change of HQs momentum from 
p to p — k is given by, 

/d 3 q ^ 

J^yf'^ va P ,q^p—k,q+k (3) 

where /' is the phase space distribution of the particles in the bulk, v is the relative 
velocity between the two collision partners, a denotes the cross section and Qq is the 
statistical degeneracy of the particles in the bulk. 

Considering only the soft scattering approximation [37], the integro-differential 
Eq. 0 can be reduced to the Fokker Planck equation: 


A(p)f 


dl = d_ 

dt dpi 
where the kernels are defined as 




and 


Ai= d 3 kw(p,k)ki 


Bij = — J d 3 kw(p,k)kikj 


(4) 


(5) 


( 6 ) 


In the limit | p |— y 0, Ai —y yp* and B t j — y DSij where 7 and D are the drag and diffusion 
coefficients respectively. We notice that under the assumption of soft collisions, the 
non-linear integro-differential equation, Eq. |T] reduces to a much simpler linear partial 
differential equation, Eq. [4] provided the function f'(q) is known. 

Now we evaluate the 7 and D for HQs interacting elastically with the bulk particles 
in the pre-equilibrium phase that appear in HIC before thermalization. The 7 for the 
process HQ[pf) + g(p 2 ) HQ(p 3 ) + g(pf) (pf s are the respective momenta of the 

colliding particles) can be expressed in terms of Ai [3?J (see also p 8 ] [39] 00]) as: 

7 = piAi/p 2 , 


(7) 
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where Aj is given by 

= 1 r d 3 p 2 r d 3 p 3 r d 3 p A 

^ 1 ~ 2E P1 J (2 t t) 3 E P2 J (2n) 3 E P3 J (27 t) 3 E P4 

— '52\M\ 2 (2Tr) 4 8 /L (p l +p2-p 3 ^,pd 

9hq 

/(P2){1 ± /(P4)}[(Pl - Ps)i] = (((Pi - Ps)i))- 

( 8 ) 


where guQ denotes the statistical degeneracy of HQ, /(p 2 ) is the momentum 
distributions of the incident particles, 1 ± fipf) is the momentum distribution (with 
Bose enhancement or Pauli suppression) of the final state particles in the bath and 
\M\ 2 represents the square modulus of the spin averaged invariant amplitude for the 
elastic process, evaluated here using pQCD. The drag coefficient in Eq. (jHJ) is the measure 
of the average of the momentum transfer, p\ — p$, weighted by the interaction through 
|A4| 2 . Similarly the momentum diffusion coefficient, D, can be defined as: 


D = - 
4 


«P3» 


«(ffi-P3) 2 » 


Pi 


(9) 


and it represents the averaged square momentum transfer (variance) through the 
interaction process mentioned above. 

The following general expression has been used to evaluate the drag and diffusion 
coefficients numerically with an appropriate choice of Z{p), 


•C Z(p) »= 


1 1 r°° p 2 dp 2 d(cosx) 

5127r 4 E Pl Jo E P2 


/(p 2 ){l ± f(P4)}'- 


A 2 (s,m 2 m 2 ) /■-! 




d(cos9 c . r 


1 

9hq 


_ r 2n 

J2\M\ 2 d<p c . m Z(p) 
Jo 


( 10 ) 


where = {s - (% +m M ) 2 }{s- (m pt -m pa ) 2 }. 


3. Initial conditions 

In most of the earlier works the distribution functions in Eq. (fill are taken as equilibrium 
distribution, i.e. Fermi-Dirac for quarks and anti-quarks and Bose-Einstein for gluons. 
The transport coefficients are then corresponding to the motion of the HQs in a 
thermalised medium assumed to be formed within the time scale ~ 1 fm/c. 

In this article instead, we will evaluate the drag and diffusion coefficients of HQ 
propagating through a non-thermal gluonic system. We achieve this by substituting 
the distribution functions in Eq. (ITU]) by pre-equilibrium distribution of gluons to be 
specified later. We choose the normalization of the non-equilibrium distributions in such 
a way that the gluon density for the chosen distribution and the thermal distribution 
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at initial temperature coincides. The initial temperature, T t is taken as 0.34 GeV for 
RHIC and 0.51 GeV for LHC collision conditions. These initial temperatures are chosen 
from simulations [Ml 1351 |36] done to reproduce V 2 and spectra of light hadrons, as 
well as the Raa and u 2 of heavy mesons [12, 33]. We compare the drag and diffusion 
coefficients of HQs in the pre-equilibrium era with those in the thermalized era at the 
initial temperatures just mentioned. 

We now specify the out-of-equilibrium gluon distribution used in this work for 
evaluating HQs drag and diffusion coefficients. According to the general understanding 
of the dynamics of pre-equilibrium stage the initial strong gluon fields (the glasma) 
shatters into gluon quanta in a time scale which is of the order of the inverse of the 
saturation scale Q s \ therefore any model of the pre-equilibrium stage will include a Q s . 
The first one we consider is the classical Yang-Mills (CYM) gluon spectrum (see [4TJ 
for details), which assumes the initial gluon fields can be expanded in terms of massless 
gluonic excitations. 

Beside CYM we also consider the model known as factorized KLN model [121113] 
which includes the saturation scale in an effective way through the unintegrated gluon 
distribution functions which, for the nucleus A, that participates in A + B collision 
reads: 



( 11 ) 


a similar equation holds for nucleus B. The momentum space gluon distribution is then 
given by 



( 12 ) 


where A f is an overall constant which is fixed by the multiplicity. 

4. Results and discussions 

We have evaluated the drag and diffusion coefficients of HQs propagating through 
a system of out-of-equilibrium gluons formed at the very early era of HIC at the 
RHIC and LHC energies. The magnitudes of the transport coefficients evaluated in 
the pre-equilibrium era are compared with those obtained in the equilibrium phase 
(both kinetic and chemical) keeping the number of particles fixed in both the cases as 
mentioned above. The temperature dependence given in Ref. |44j has been used to 
estimate the value of the strong coupling, a s at T = T t . The Debye screening mass, 
mD = ^8a s (N c + Nf)T 2 / tt, for a system at temperature, T with N c colors and Nf 
flavors has been used to shield the infra-red divergence associated with the f-channel 
scattering amplitude, mo is estimated at T = 0.34 GeV for RHIC energy and at 
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Figure 1 . Drag coefficient (left panel) and diffusion coefficient (right panel) as a 
function of momentum for charm quark at RHIC energy. The results corresponding 
to thermal gluons (kinetic equilibrium) and thermal quarks and gluons (chemical 
equilibrium) are evaluated at a temperature 340 MeV. 


T = 0.51 GeV for LHC energy. For the pre-equilibrium and equilibrium system same 
values of mr_> and a s have been used. Later in this section we will show results where 
mu is estimated in a self-consistent way with the underlying distribution function. 

In the left panel of Fig. CD we plot the drag coefficient of the charm quark as a function 
of momentum in the pre-equilibrium phase for CYM and KLN gluon distributions and 
compared the results with the equilibrated phase (both kinetic and chemical) at T=0.34 
GeV. We find that the magnitude of the drag coefficient in the pre-equilibrium phase is 
quite large, indeed comparable to the one in the kinetic equilibrium phase, indicating 
substantial amount of interaction of the charm quarks with the pre-equilibrated gluons. 
We notice that 7 for the case of chemically equilibrated QGP (dotted line) is smaller 
than the one of the purely gluonic system (solid line in Fig. [[]). This can be understood 
by considering the fact that to keep the total number of bath particles fixed some of 
the gluons in the purely gluonic system has to be replaced by quarks in the chemically 
equilibrated QGP. Because the gluon appears in more colours than quarks the cross 
section for eg interaction is larger than cq which leads to larger drag of HQs in a 
pre-equilibrated gluonic system than in a chemically equilibrated QGP. The CYM 
distribution gives larger value of the drag coefficient compared to KLN distribution. The 
KLN provides harder momentum distribution (compared to CYM) hence have a smaller 
difference with the HQs distribution. Therefore, the momentum transfer between the 
gluon (with KLN distribution) and the HQ is small. As the drag coefficient is a measure 
of the momentum transfer weighted by the interaction strength, the KLN gives rise to 
lower drag compared to the one obtained with the CYM distribution. 

The variation of D with momentum for charm quarks has been depicted in the 
right panel of Fig. [D We find that the magnitude of the D for CYM initial condition is 
similar to the case where the gluons are kinetically equilibrated. However, the D for the 
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Figure 2. Dj 7 as a function of momentum for charm quark at RHIC energy. The 
result corresponding to equilibrium cases are evaluated at a temperature 340 MeV. 


KLN initial condition is larger than CYM. As mentioned before, the KLN has harder 
momentum distributions than the CYM, hence KLN distribution has larger variance 
compared to CYM. Since the momentum diffusion is a measure of the variance acquired 
through interaction with the bulk, the charm diffusion coefficient is larger for KLN 
distribution, which is reflected in the results displayed in Fig. [0 

In Figj2] we have depicted the diffusion to drag ratio, D/ 7 as function of 
momentum. D /y can be used to understand the deviation of the calculated values 
from the value obtained by using Fluctuation-Dissipation theorem (FDT) (green line in 
the figure). Since the KLN has a harder momentum distribution, results obtained from 
KLN input deviates from FDT more. 

In this work, the dependence of the drag/diffusion coefficient on the variation of the 
phase space distribution is addressed. To make the present study more consistent the 
effects of phase space on the dynamics through Debye screening mass has been included. 
It is well-known that the Debye mass depends on the equilibrium distribution as follows: 

m 2 D = Tra s g G j J^h N cfg + N ff q ) ( 13 ) 

It is interesting to see how the results are affected when the equilibrium distribution in 
Eq. [13] is replaced by the KLN or CYM distributions. 

Replacing the thermal distribution of gluons in Eq. (fT3j) by the KLN and CYM 
gluon distributions and setting f q — 0 for a gluonic system we estimate m D and use 
this to calculate 7 and D. The momentum variation of 7 and D are displayed in Fig. [3] 
for charm quarks. It is interesting to note that in this case 7 is almost unaffected by 
the bulk distributions, i.e. of KLN, CYM and thermal gluons. We End some difference 
between the aforementioned bulks and the chemically equilibrated QGP, which is caused 
by the less number of gluons in the latter case as discussed earlier. 

The drag (left panel) and diffusion (right panel) coefficients for b quarks in the 
pre-equilibrium phase have been displayed in Fig. [4] We notice that these coefficients 
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Figure 3. Drag coefficient (left panel) and diffusion coefficient (right panel) as a 
function of momentum for charm quark at RHIC energy. Debye mass is computed 
self consistently using Eq. [T3] The results corresponding to thermal gluons (kinetic 
equilibrium) and thermal quarks and gluons (chemical equilibrium) are evaluated at a 
temperature 340 MeV. 




Figure 4. Drag coefficient (left panel) and diffusion coefficient (right panel) as a 
function of momentum for bottom quark at RHIC energy. Debye mass is computed 
selfconsistently using Eq. [13] The results corresponding to thermal gluons (kinetic 
equilibrium) and thermal quarks and gluons (chemical equilibrium) are evaluated at a 
temperature 340 MeV. 


are smaller for b than c. However, the qualitative variation of b diffusion coefficient with 
momentum is similar to c. 

In the left panel of Fig. [5] the momentum variation of the 7 of the c quark in the 
pre-equilibrium phase is depicted for LHC collision conditions. This result is compared 
with the one computed for a QGP at temperature T=0.51 GeV. Here same values of the 
coupling for both the equilibrium and pre-equilibrium system have been used. The value 
of mp is taken from Eq. dT3]l for all the cases. We find that the magnitude of the drag 
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Figure 5. Drag coefficient (left panel) and diffusion coefficient (right panel) as a 
function of momentum for charm quark at LHC energy. Debye mass is computed 
selfconsistently using Eq. [Hi] The results corresponding to thermal gluons (kinetic 
equilibrium) and thermal quarks and gluons (chemical equilibrium) are evaluated at a 
temperature 510 MeV. 




Figure 6. Drag coefficient (left panel) and diffusion coefficient (right panel) as a 
function of momentum for bottom quark at LHC energy. Debye mass is computed 
selfconsistently using Eq. [13] The results corresponding to thermal gluons (kinetic 
equilibrium) and thermal quarks and gluons (chemical equilibrium) are evaluated at a 
temperature 510 MeV. 


coefficient in pre-equilibrium phase is similar to that obtained for the kinetic equilibrium 
case, however, for a QGP the value of drag is smaller because of the less number of gluons 
as mentioned earlier. In the right panel of Fig. [5] we show the variation of D with p 
for the c quarks at LHC energy. Again the variation of the diffusion coefficient at LHC 
energy is qualitatively similar to RHIC. Similar to RHIC the diffusion at LHC conditions 
is larger for the KLN distribution. 

Similarly the 7 and D for b quarks are plotted in Fig. [ 6 ] for LHC collision conditions. 
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Qualitatively the variation of these coefficients with p is similar to RHIC. Although, the 
quantitative values at LHC collision conditions are larger than RHIC as expected owing 
to the larger density and temperature.. 

5. Summary and discussions 

In this work we have studied the momentum variation of the drag and diffusion 
coefficients of HQs in the pre-equilibrium era of heavy ion collisions at RHIC and LHC 
energies. The momentum distribution of the pre-equilibrated gluons have been taken 
from CYM and KLN formalisms. This study is motivated by the fact that the effect of 
the pre-equilibrium phase on the HQs suppression and elliptic flow might be relevant 
for low-energy nuclear collisions. For example, the simulations of Au+Au collision at 
RHIC energy show that the equilibration is achieved approximately within a time scale 
of 1 fm/c, while the lifetime of the QGP phase turns out to be about 5 fm/c [HU 135] : 
hence the system spends about 20% of QGP life-time in the pre-equilibrium phase. In 
the case of Pb-Pb collisions at LHC energy the equilibration time is shorter and lifetime 
of the QGP phase is larger, hence in this case we expect a smaller effect vis-a-vis pre¬ 
equilibrium phase. 

We have compared the magnitudes of the transport coefficients computed for 
equilibrated and pre-equilibrated system, keeping the number of particles same in 
the two cases. We have found that the magnitude of the transport coefficients in 
the pre-equilibrium phase is comparable to the values obtained with a thermalized 
gluonic system. Moreover, we have also found that the transport coefficients in the 
pre-equilibrium era are larger than the ones obtained for a chemically equilibrated QGP 
system. This is due to the fact that for a fixed number of particles the number of 
gluons are less in the equilibrated QGP than the pre-equilibrated gluonic system and 
the HQ + q cross section is smaller than the HQ + g cross section. 

The results obtained in this work may have significant impact on the experimental 
observables, for example on heavy mesons Raa and elliptic flow, as well as on the 
suppression of single electron spectra originating from the decays of heavy mesons and 
their elliptic flow. We will address these aspects in future works. 
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